Persistent currents on graphs 
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We develop a method to calculate the persistent currents and their spatial distribution (and trans- 
port properties) on graphs made of quasi-lD diffusive wires. They are directly related to the field 
derivatives of the determinant of a matrix which describes the topology of the graph. In certain 
limits, they are obtained by simple counting of the nodes and their connectivity. We relate the 
average current of a disordered graph with interactions and the non-interacting current of the same 
graph with clean ID wires. A similar relation exists for orbital magnetism in general. 

PACS Numbers: 73.35, 72.20M, 73.40L 
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The existence of persistent currents in mesoscopic 
metallic rings is a thermodynamic signature of phase 
coherence Q. These currents have been calculated us- 
ing diagrammatic methods in which disorder and inter- 
actions are treated perturbatively in a way very 
similar to the calculation of transport quantities like the 
weak-localization (WL) correction, or the universal con- 
ductance fluctuations (UCF). Like transport quantities 
, they have also been derived (after disorder averag- 
ing) using semi-classical calculations, in which they were 
expressed in terms of the classical and interference parts 
of the return probability for a diffusive particle [p|-p^[. 
This formalism had made possible the calculation of WL 
corrections on any type of graph made of diffusive wires 
p3[ . Diffusion equation was solved on each link of the 
graph with current conservation on each node. For a 
network with N nodes, the return probability could be 
related to the elements of a N x N " connectivity" matrix 
M and its inverse. This method has also been used re- 
cently to calculate the magnetization of such a network 
Ej l, but it required rather lenghty calculations. 

In this letter, we show that the magnetization and the 
transport quantities can be directly written in term of 
the determinant det M of the connectivity matrix. Be- 
sides being a very powerful method to calculate the above 
quantities, this result leads to a straightforward harmonic 
expansion of these quantities for any network geometry. 
The efficiency of this method is shown for simple geome- 
tries of connected rings. In addition, we are able to derive 
the local distribution of the currents in the links of the 
network. Since the persistent current problem has still 
to be considered as unsolved, it is of interest to motivate 
new experiments in various geometries for which the mag- 
netization and its distribution can be simply predicted 
and related to geometrical or topological parameters. 

In the course of this work , we shall obtain a simple 
expression for the spectral determinant of the diffusion 
equation, defined as: 



where E n are the eigenvalues of the diffusion equation 
and b n are regularization factors Jl5| . Using the analogy 
between the diffusion and the Schrodinger equation, we 
will point out a very simple relation between the Hartree- 
Fock (HF) average magnetization of a diffusive system 
and the grand canonical magnetization of the correspond- 
ing clean system. As a simple example, we relate the 
Aslamasov-Larkin contribution to the magnetization and 
the Landau susceptibility. 

All quantities of interest in this work can be related 
to the solution P(f,r' ,uj) of the diffusion equation in a 
magnetic field B — V x A(f) = 1 throughout the 

paper): 

iu + j - D(Vp-2ieA) 2 ] P(r,r',uj) =5(f f -r') (2) 

D is the diffusion constant. Unless specified, the mag- 
netic field dependence is implicit. 7 = l/r^ = D/L 2 ^ is 
the phase coherence rate. and are respectively 
the phase coherence length and time. In the follow- 
ing, we will only need the space integrated return prob- 
ability P(t) = f d d rP(r,r,t). It is simply written in 
terms of the eigenvalues E n of the diffusion equation, 
p ( f ) = E„ e ~ (B " +7) * = fb(*)e~ 7 *. The time integral of 
P(t), i.e. the Laplace transform of Po(t) can be straight- 
forwardly written in terms of the spectral determinant 
(0). 
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Let us now recall how average magnetizations can be 
written in terms of P(t). Here we restrict ourselves to 
T = OK. The fluctuation of the magnetization Mt yp = 
((M 2 ) - (A/) 2 ) 1 / 2 are given by 0: 



+ - dt ^m^m (4) 



Ml) =Y[b n {j + E n ) 



(1) 



where P"(t,B) = d 2 P(t, B)/dB 2 . The main contribu- 
tion to the average magnetization is due to electron- 
electron interactions ||,|J . Considering a screened inter- 
action U(r — r') — US(f— f') and defining Ao = Upo 
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where po is the density of states (DoS) at the Fermi en- 
ergy cf, the Hartree-Fock (HF) contribution to the mag- 
netization has be written as [Gol: 



(M ee ) = - 
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Considering higher corrections in the Cooper channel 
leads to a ladder summation [[l8 17 5|J[^], so that Ao 
should be replaced by X(t) = A /(l + A ln(e F i)) [^9). 
We shall discuss later the contribution of this renormal- 
ization. 

Using standard properties of Laplace transforms, the 
above time integrals can be written as integrals of the 
spectral determinant, so that the magnetizations read: 
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In the case of a ring or a graph geometry, the integral 
converges at the upper limit. For the case of a magnetic 
field in a bulk system, this limit should be taken as l/r e 
where r e is the elastic time. Finally, we also recall that 
transport properties like WL or UCF can be also related 
to the spectral determinant jL4j. 

We now wish to emphasize an interesting correspon- 
dence between the HF magnetization of a phase coher- 
ent interacting diffusive system and the grand canonical 
magnetization Mq of the corresponding non-interacting 
clean system. The latter can also be written in term of 
a spectral determinant. The grand canonical magnetiza- 
tion Mq is given quite generally by: 



M 



= / deN(e) 



dB 



dB 



(8) 



where the integrated DoS is 



N(e) = --Im Vln(e M - e+) = --ImlnS(e + ) (9) 

TT ' "7T 



where e + = e + iO, S(e) = JT e(I M e M ~ e ) = 5rf (7 = ~ e )- 
where are the eigenvalues of the Schrodinger equation. 

For a clean system these eigenvalues are the same as 
those of the diffusion equation, with the substitutions 
D -> h/(2m) and 2e e @. 

Comparing eqs.(||,||) with eq.(0), we can now formally 
relate Mo and the HF magnetization (M ee ) of the same 
diffusive system: 

M = -lim Ao ^OT^Im[(M ee )(-e F - jo)] (10) 
Ao 

This limit corresponds to taking the first order contribu- 
tion in Ao- As a simple illustration, consider the orbital 
magnetic susceptibility of an infinite disordered plane. 



For a disordered conductor, it is the Aslamasov-Larkin 
susceptibility \al Jig]: 
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To = e F e 1 / A ° and cf>o = h/e is the flux quantum. After 
replacing 7 by — ep — *0, taking the imaginary part of the 
logarithm and replacing D and 2e, we recover the Landau 
susceptibility for the clean system: \o = — e 2 /(247rm). 

We now calculate the spectral determinant for quasi- 
1D graphs. By solving the diffusion equation on each 
link, and then imposing Kirchoff type conditions on the 
nodes of the graph, the problem is reduced to the solution 
of a system of N linear equations relating the eigenvalues 
at the N nodes. Let us introduce the N x N matrix M 
& 



M aa = 2jcoth(?? a ff) , M a p = 

P 



sinh r] a t 



(12) 



The sum extends to all the nodes (3 connected to 
the node a; lap is the length of the link between a and 
P- T) a p — l a p/L,p. The off-diagonal coefficient M a p is 
non zero only if there is a link connecting the nodes a 
and /3. 9 a p = (47r/0o) j A.dl is the circulation of the 
vector potential between a and f3. The authors of ref. 
|l3| derived a relation between V and the elements of the 
matrix M and its inverse T = M _1 : 



2 1 V = (N -N B ) + J2 V*pF a p 

(a/3) 



F aP = cothry Q/3 - ^+ T ^ + 2Re(e w ^Tp aJ 



(13) 
cosh 770,^3 



sinh rj af: 



sinh rj a 







where Nb is the number of links in the graph. Using the 
equality: Tr(M _1 <9 7 M) = <9 7 lndet M and recognizing in 
each term of (^) the partial derivative with respect to 
7, we find that eq.(|l3|) can be rewritten as: V — -g- \nSd 
where the spectral determinant Sd is given by: 



S d 



L 



N B -N 

) TT sinh?7 Q/ 3 detM 



(14) 



apart from a multiplicative factor independent of 7 (or 
L^). Lq is an arbitrary length. We have thus transformed 
the spectral determinant which is an infinite product in 
a finite product related to det M. 

As an example, we consider a disordered ring of 
perimeter L, to which one arm of length b is attached. 
The spectral determinant is equal to: 

Sd = sinh Ry sinhy + 2coshi?y (coshy — cos(47r</?)) 
where ip = 4>/4>o is the ratio between the flux <fi threading 
the ring and the flux quantum, y = (L/ L^) and R = b/ L. 
Thus the average magnetization is: 
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(M ee ) = ^ 



2 sin 4-7T(/? ydy 



tanh i?y sinh y + 2 (cosh y — cos 4inp) 

(15) 



If there is no arm (R = 0), we retrieve the classical 
expression for the average magnetization of a disordered 
ring p2]j . We notice that, in the limit L^, the mag- 
netization remains finite and is equal to 2/3 of the single 
ring magnetization (for < L, which corresponds to 
typical experimental values). 

We want first to outline once more the connection be- 
tween ballistic and disordered regimes. From eq. @ and 
with the mapping (|To|), 7 — > —£7 — iO and L/L^ — * 
ikL where k is the wave vector of the solutions of the 
Schrodinger equation, we immediately recovers the cur- 
rent in a one channel ballistic ring p3[ . 

Let us come back to a diffusive network made of con- 
nected rings. Experimentally, the coherence length is of 
the order of the perimeter of one ring so that only a few 
harmonics of the flux dependence may be observed. It is 
then useful to make a perturbative expansion. We split 
the matrix as M — D — N, where D is a diagonal matrix: 
D aa — Maa ~ z a to the lowest order in (z a is the con- 
nectivity of the node a); N a/3 = M a p ~ 2e _i °' ,/I "*e if, =< 3 . 
Expanding bidet (J - D~ l N) = Tr[ln(J - D^N)], we 
have: 



IndetM = lndetD- ^ -Tr[(Zr 1 .ZV) n ] 



(16) 



We call "loop" I, a set of n nodes linked by n wires in a 
closed loop. The length Li of a loop I is the sum of the 
lengths of the n(l) links. The flux dependent part of In S 
can be expanded as: 



. e -Li/L* cos (4 7r i /0 o ) (17) 



is the flux enclosed by the loop /. 
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FIG. 1. connectivity factors (2/ zi){2/ Z2) entering in the 

loop expansion (|l7|), for a series of identical connected rings, 
a single ring, and a ring with one arm. 

For example, we consider the cases shown on Fig j^. 
Reducing the above sum to elementary loops Iq (with 
two nodes), so that n(lo) = 2, the first harmonics of the 
total magnetization, to the first order in Ao is: 

(M ee ) = 2G^(L/L + l)e- L ' L * (18) 



over the m rings of the structure (see Fig.[j]). In par- 
ticular, it is G — (m + 2)/4 for an open necklace of m 
rings and G = m/4 for a closed necklace. The same re- 
duction factors were obtained for weak-localization cor- 
rections after lengthy calculations for m = 1, 2, 3, 00 in 
rcf. p3| . For the isolated ring, one recovers the known 
first harmonics |ll]] and the above reduction factor 2/3 
for the ring with one arm. For an harmonic p of the 
magnetization, corresponding to a winding number p in 
the diffusion process, one should renormalize the inter- 
action parameter because of the Cooper renormalization 
A = Ao/(l + A lne F /(7r£ c /p 2 )) [§. 

Fig.|2j displays a comparison between the magnetization 
of different networks of connected rings, evaluated nu- 
merically using eqs. (|?],|l4|) . The perturbative expansions 
are in extremely good agreement with exact results as 
soon as the coherence length is smaller than the perime- 
ter of one ring (see dashed lines in FigJ|). 
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FIG. 2. Magnetization per ring for networks of connected 
rings normalized to the single ring magnetization, calculated 
exactly (solid lines) and with the loop expansion (dashed 
lines). The perimeter of all rings, and side-arm lengths are 
equal to L, The three bottom curves correspond to regular 
networks made of an infinite number of rings (only 3 are repre- 
sented). In these cases, the magnetization has been divided by 
the number of rings. The flux threading all rings is (f> — </>o/8 

Finally, we calculate the distribution of the local cur- 
rent on each link of the graph. On a link {a(3) , the aver- 
age current is given by the derivative of the Hartree-Fock 
energy correction Ejjf to the vector potential A(r), where 
r is any point belonging to the link (a/3): 
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(19) 



J^M) = ^MM^M a0 ) (20) 



Fluctuations of the current corresponding to eq. 
be obtained similarly p4| . 



can 



where G = X]{z } ^/( z i z 2)- %\ and Z2 are the connec- 
tivity of the two nodes of each loop. The sum is made 
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FIG. 3. Current distribution (left: average current; right: 
variance) for a square network, in the limit < L, i.e. when 
the flux dependence of the current is harmonic. The numbers 
show the amplitude of the average and typical magnetization 
per plaquette, in units of the magnetization of the single ring. 
It is maximum at the corner plaquettes. The thickness of each 
link is proportional to the amplitude of the current on this 
link, obtained by difference (resp. sum) of the average (resp. 
typical squared) magnetizations of the plaquettes neighboring 
the link. 

In the limit < L considered above, the current dis- 
tribution can also be derived quite simply. Indeed, in this 
approximation, the total magnetization can be written as 
a sum (M ee ) ~ X)fc( TO fe) where (rrik} is identified as the 
magnetization of a plaquette k and depends on the posi- 
tion of this plaquette in the array. It is given by the rules 
of eq. (O) and is shown on Fig|| for a regular square lat- 
tice. The average persistent current flowing in one link 
is the difference of the two plaquette currents neighbor- 
ing it. The distribution of average current is sketched on 
Fig§ 

The fluctuations can be described in the same way: 
namely it is a sum of terms which can be interpreted as 
fluctuations of the magnetization of one plaquette: thus 
the fluctuations of plaquettes are independent, and the 
fluctuations of current in one link is the sum of the fluc- 
tuations of its two nearby plaquette-current. 

In conclusion, we have developed a formalism which 
relates directly the persistent current, and the transport 
properties (although not detailed in this letter) to the de- 
terminant of a matrix which describes the connectivity 
of the graph. From a loop expansion of this determinant, 
simple predictions for the magnetization and the spatial 
distribution of the persistent current in any geometry can 
now be compared with forthcoming experiments on con- 
nected and disconnected rings. We have also found a 
correspondence between the phase coherent contribution 
to the orbital magnetism of a disordered interacting sys- 
tem and the orbital response of the corresponding clean 
non-interacting system. 

We acknowledge useful discussions with E. Akkermans, 
A. Benoit, E. Bogomolny, H. Bouchiat, D. Mailly and L. 
Saminadayar. 
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